Qubit measurements with a double-dot detector 
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We propose to monitor a qubit with a double-dot (DD) resonant-tunneling detector, which can 
operate at higher temperatures than a single-dot detector. In order to assess the effectiveness of 
this device, we derive rate equations for the density matrix of the entire system. We show that the 
signal-to-noise ratio can be greatly improved by a proper choice of the parameters and location of 
the detector. We demonstrate that quantum interference effects within the DD detector play an 
important role in the measurement. Surprisingly, these effects produce a systematic measurement 
error, even when the entire system is in a stationary state. 
PACS: 73.50.-h, 73.23.-b, 03.65.Yz. 



The single electron transistor (SET) is a sensitive de- 
vice for quantum measurements [1-3]. It can be used as 
a monitor of a charge qubit, provided that the energy 
level Eq carrying the current is close to the Fermi levels 
of the reservoirs ^j-l.r [see Fig. 1(A,A')]. Then due to the 
electrostatic repulsion U between the electrons, the SET 
current drops when the qubit is in the state E2, as in 
Fig. l(A'). 

It is clear that one needs very low reservoir temper- 
atures in order to use the SET as a sensitive detector. 
This requirement, however, can be weakened by tak- 
ing a double-dot (DD) for monitoring the qubit state. 
Fig. 1(B,B'). Similar to the SET, the DD current de- 
creases sharply whenever the electron of the qubit is close 
to one of the dots, Fig. l(B'). In contrast with the SET, 
however, the reservoirs' temperature (T) will not affect 
the current li - T :$> Eq :$> ^lr + T [4,5]. 
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Fig. 1: The qubit measurements using the SET de- 
tector (A,A') and the DD detector (B,B'). The detector 
current / drops down whenever the state E2 of the qubit 
is occupied. F^ represent tunneling rates to the reser- 
voirs and 7 denotes the interdot coupling. 

Because of quantum interference effects, the dynam- 
ics of the measurement process using the DD detector is 
more complicated than with the SET detector. An elec- 



tron flowing through the DD can be trapped in a linear 
superposition of the dot states. As a result, quantum in- 
terference could modify the signal in such a way that the 
DD cannot monitor the qubit. It is necessary to analyze 
the influence of the DD on the qubit motion (and vice 
versa) in order to establish the optimal conditions for uti- 
lizing the DD as an effective quantum detector. This can 
be done by solving the Schrodinger equation describing 
the combined system of qubit and detector. 

In fact, the setup shown in Fig. 1(B,B') represents a 
generic class of non-demolition quantum measurements 
where a measured system interacts with only one state 
of the apparatus, while the apparatus may be in a super- 
position of states. This can take place in many devices 
based on interference, for instance in Electronic Mach- 
Zehnder Interferometer [8,9]. We therefore believe that 
our analysis of a qubit interacting with the DD detector 
can be useful for understanding many different quantum 
measurements. 

Let us describe the entire setup shown in Fig. 1(B,B') 
by the tunneling Hamiltonian H = Hq -\- Hdd + -ffmt, 
where 

Hq = Eia\ai + E2a\a2 + ^{a\a2 + a\ai) , 
Hdd = Ho + Eo{c\ci + c\c2) + 7(clc2 + c\cx) 
\iSAcXc\ + r!f 4cf + H.c.) + C/i24ci4c2 



Hint = J7a2a24c2 



(1) 



are the qubit and the DD Hamiltonians, and Hint is 
their interaction. Here a}(a) is the creation (annihila- 
tion) operator for the electron in the qubit and c^(c) 
is the same operator for the DD; is the coupling be- 
tween the states |a| 2IO) of the qubit, and 7 is the cou- 
pling between the states |c| 2IO) of the DD. The Hamil- 
tonian Hq = Y.x[HxiciVc{ + E^ic^Yc^] describes the 
reservoirs, where fl^' are the couplings between the 
right and left dots with the right and left reservoirs. 
We assume weak energy dependence of these couplings, 

L R 

fl^' ~ flL,R- Then the corresponding tunneling rates 
are Tl^h = 2TTpL^i^flj^ ^, where pl,r are the density of 
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states in the reservoirs. The latter quantities arc also 
weakly dependent of energy. The last term in de- 
scribes the inter-dot repulsion. For simplicity we consider 
electrons as spinlcss fermions. 

Using the technique developed in Ref. [10,11] for the 
case of large bias voltage, V = — Hr, we can partially 
trace out the reservoir states in the equation of motion 
for the density matrix of an entire system, ig = [H,q]. 
As a result we arrive at the Bloch-type rate equations for 
the reduced density-matrix, afj{t), describing the qubit- 
detector evolution, where the indices i,j denote all avail- 
able discrete states of the detector-qubit system and n 
denotes the number of electrons which have arrived at 
the right reservoir by time t. In our case, Fig. 1(B,B'), 
the available discrete states are labeled (a, b, c, d), denot- 
ing the cases that the DD is empty (a), the left dot of 
the DD system is occupied (&), the right dot of the DD 
system is occupied (c), and both dots are occupied (d), 
while the electron of the qubit occupies the level Ei (see 
Fig. IB). Correspondingly, {a' ,b' , c' , d') denote the same 
states but where the electron of the qubit occupies the 
level E2 (see Fig. IB'). If the inter-dot repulsion is large, 
U12 ^ V, the states d, d' do not contribute terms to the 
equations of motion. We obtain in this case [10-12]: 
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Note that these equations are obtained from the origi- 
nal many-body equations ig = [H, g] without the explicit 
use of any Markov-type or weak-coupling approximations 
in the case of large bias voltage, V » ^l,r, U [10,11]. 
There are no other limitation on U, in contrast with our 
analysis of the SET detector [12]. 

Equations (2) are different from the standard master 
equations, describing quantum system interacting with 



the environment (detector) by keeping track of the en- 
vironment variables. In our case this is the number 
of electrons (n) arriving the collector. This allows us 
to find the time evolution of the qubit and the detec- 
tor at once. For instance, the qubit behavior is de- 
scribed by the (reduced) density matrix aq{t) = {<Ta0{t)} 
with = {1,2}, where an = EJ<a + ^Ib + <c), 

C^12 = En«o' + <6' + <c') and (722 = 1 - 0-11. 

On the other hand, by tracing out the qubit variables 
wc obtain the probability of finding n electrons which 
have arrived at the collector, P„(i) = J2j'^jj{i)- This 
quantity allows us to determine the average detector cur- 
rent and its shot-noise spectrum. The former is given by 



I{t) = e^nP„(t) = eTRaRit), 



(3) 



where aR{t) = J2nWcc{t) + '^"c'(*)] ^lie probabihty 
that the right dot is occupied. The shot-noise spectrum, 
S{w), is obtained from the McDonald formula [13,14] 

S{uj)=2e'^oj (itsin(a;t) Vn2p„(t), (4) 

One finds from Eqs. (2), (4) that 

S{lo) = 2e''uTRlm [Zec(w) + ^c'c' (w)] , (5) 

where Zy (w) = Y.ni'^n + l)(j^j{t) exp{iujt)dt. These 
quantities are obtained directly from Eqs. (2) which are 
reduced to a system of linear algebraic equations after 
the corresponding integration over t [15]. 

Consider first the static qubit, f2 = 0. Solving Eqs. (2) 
for this case one finds that the stationary current, / = 
I{t 00) obtains the value Ti = TRaR{U = 0) when the 
qubit is in the state Ei, and I2 = ^rO'r{U) when the 
qubit is in the state E2, Fig. 1(B,B'), where 



^r{U) 



(6) 



As expected, the detector current decreases whenever the 
electron of the qubit is close to the DD detector. Con- 
sider now 7^ 0. Wc assume that for an "ideal" detector 
its average current would follow the qubit motion [16] , 



I{t)=haxi{t)+l2[l-<Jri{t)]. 



(7) 



This condition, however, cannot be fully met since the de- 
tector's response is limited by the rate of tunneling from 
the right dot to the collector. Nevertheless, if this tran- 
sition is fast enough compared to the qubit frequency, 
Ffl J> O, one expects to approach Eq. (7) . 

Let us compare <Tii(t) with the average "signal," [/(t) — 
J2]/AJ, where A/ = Ii — I2. The results of our calcu- 
lations for 7 = O are presented in Fig. 2. The initial 
conditions correspond to the qubit electron in the upper 
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dot and the detector current I{t = 0) = h- One finds 
that the detector does not follow the qubit oscillations 
well when Tu = 7, Fig. 2(a). On the other hand, in 
Fig. 2(b) where Tu 7, the detector performance is 
much improved, in accordance with our arguments. 





(a) (b) 

Fig. 2: The probability of finding the qubit in the state 
El, (dashed line) compared with the average detector 
signal, [I{t) — /2]/A/ (solid line) as a function of time 
for 7 = f7 and U = 5f7, = 5f7. 

The results displayed in Fig. 2(a) are surprising. One 
expects that in the steady-state limit (t ^ 00) the aver- 
age detector current should be distributed between the 
values Ii and I2, with probabilities an and I — an to find 
the qubit in the states -E1.2, respectively. Equation (7) 
should thus always hold in the limit of i — > 00 for any 
such device (for instance, the SET detector [12]). In the 
case of the DD detector, however, Eq. (7) does not hold 
in the steady-state limit, as seen in Fig. 2(a). In fact, this 
can be obtained analytically in the limit of small U by 
expanding the stationary current, / = I(t ~> 00), Eq. (3), 
in powers of U. One finds for the detector's signal: 
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It follows from this expression that a mismatch between 
the signal and the qubit (dn) survives even in the limit 
U ^ 0. In this case it depends only on the ratio Tj^/il. 
The other detector parameters 7 and F^ enter only in 
the term proportional to U^. 

So where is the "hidden" probability that is responsi- 
ble for the systematic error in the qubit measurements? 
It can be recovered in the linear superposition of the de- 
tector and qubit states. The DD current flows via two 
discrete energy levels, Eq and Eq + U. A carrier wave 
function thus proceeds through a linear superposition of 
these states, b{b') and c(c'). The qubit is itself a two 
level system described by superposition. These different 
superpositions involve the same states {b,b',c,c') of the 
entire system, and hence are entangled. This is reflected 
in the off-diagonal terms abc' and a^'c, Eqs. (2). As a 
result the superposition of qubit states (qubit 's "phase") 
directly affects the DD dynamics, leading to a violation 
of Eq. (7). This would happen even in the limit U — > 00. 
In this case the state (c') disappears from Eqs. (2), but 
the off-diagonal term act' would still survive in the limit 
t — > (X). One should note that in the case of the SET such 



entanglement cannot occur. Therefore Eq. (7) holds for 
the SET, even though ai2{t ^ 00) ^ 0. 

We find from Fig. 2 and Eq. (7) that the detector's 
performance improves when Tji 3> il. At the same time, 
however, its average signal decreases. In order to as- 
sess the detector's efficiency this signal should be com- 
pared with its noise. An appropriate measure of the 
detector efficiency is the integrated signal to noise ra- 
tio [19,2], s/n = JZ,[\Is^gH\yS{Lo)]duJ/2^^. The signal 

Isig{uj) = J^[I{t) ~ I{t ^ cx))] exp(zci;i)c?i corresponds 
to a deviation of the detector current from its stationary 
value, and can be evaluated using Eqs. (2), (3), and (5). 
We show in Fig. 3 how the integrated signal to noise ra- 
tio behaves as a function of both U and the ratio Vr/^, 
for Tl — 511 and 7 = fi. A peak is observed uniformly 
throughout the range of C/ at 7 w 0.4Ffl;. The signal to 
noise ratio depends weakly on U when U/Q. > 15, allow- 
ing good operation even at low values of U . Comparing 
the performances of the DD and the SET detectors, we 
see that the maximal signal to noise ratios obtained are 
comparable. However, the SET reaches these values only 
in the asymmetric limit Vr, ^ F^ [12], while the DD 
detector's signal to noise ratio is shown to be optimized 
without such a restriction. 
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Fig. 3: The integrated signal to noise ratio of the DD 
detector. Here ^ = Q,,Tl = hVl. 

The results so far presented refer to a qubit that is 
positioned near the right dot of the detector. One may 
ask what happens when the qubit is positioned near the 
left dot. It was suggested in Refs. [20,12] that the per- 
formance of any quantum detector would improve if the 
detector were to operate mostly in the states where there 
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is no actual interaction with the qubit. By this argument 
one can expect that putting the qubit near the left dot 
would lead to poor performance. Indeed, in this case the 
misalignment of the energy levels prevents electron prop- 
agation to the right dot, so that the electron is pinned to 
the left dot. This increases the weight of states where the 
detector interacts with the qubit. On the other hand, if 
the qubit is located near the right dot as in Fig. 1(B,B'), 
the same misalignment of levels localizes the electron in 
the left dot, diminishing the occupation of the right dot, 
Eq. (6). As a result, the actual interaction with the de- 
tector decreases and so it is expected to operate better. 

Our conclusion about the asymmetry with respect to 
the qubit 's location can be confirmed by direct evaluation 
of the detector efficiency, in the same way as presented 
in Figs. 2 and 3. We can also confirm it in a different 
way by evaluating the power spectrum of the detector 
current, S{u)), via Eq. (5). This displays a pronounced 
peak at cj = 2f7, generated by the qubit oscillations. It 
was argued in Ref. [18] that the peak-to-background ra- 
tio, S{2fl)/S{u) oo), is a measure of the detector ef- 
ficiency. Figure 4 exhibits this ratio for the two qubit 
positions as a function of U. As we increase the misalign- 
ment of the levels we find that the two curves separate. 
Thus, the setup with the qubit near the right dot is in- 
deed more effective, in accordance with our arguments. 
Note that the peak-to-background ratio for this setup de- 
pends weakly on U for U/il'^ 15. We have already noted 
that the signal-to-noise ratio, which is another criterion 
of the detector efficiency, shows a similar behavior. 
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Fig. 4: The peak to background ratio of the current 
power spectrum as a function ?7 for F/, =5 0, = 10 Q 
and 7 = 0. Solid line: qubit positioned next to the 
detector's left dot, dotted: qubit positioned next to the 
detector's right dot. 

We can show that the maximal value of the peak-to- 
background ratio for the DD detector approaches 3 when 
Fij, [/ ^ (the dependence on F^ is not essential). 
The same maximal value was obtained for the SET de- 
tector [12]. Therefore, while both detectors are sensitive 
measurement devices, they do not reach the effectiveness 
of an ideal detector [18]. 

In summary, we have proposed the use of a double dot 
structure for the measurement of a charge qubit. We ob- 



tained a set of rate equations describing the entire system 
and displayed the conditions under which such a mea- 
surement is effective. We found the measurement to be 
most sensitive when the detector operates mainly in the 
states where no interaction with the qubit takes place. 
We further demonstrated that, because of quantum in- 
terference effects inside the detector, the stationary cur- 
rent is not determined solely by the probabilities of the 
stationary qubit, but reflects the qubit phase as well. 
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